arXiv:1507.01560vl [hep-ph] 6Jul2015 


Inhomogeneous phases in nonlocal chiral quark models 


J.P. Carlomagno a,b , D. Gomez Dumm a,b and N.N. Scoccola b,c,ci 
a IFLP, CONICET — Dpto. de Fisica, Universidad 
Nacional de La Plata, C.C. 67, 1900 La Plata, Argentina, 
b CONICET, Rivadavia 1917, 1033 Buenos Aires, Argentina 
c Physics Department, Comision Nacional de Energia Atomica, 
Av.Libertador 8250, 1429 Buenos Aires, Argentina 
d Universidad Favaloro, Solis 453, 1078 Buenos Aires, Argentina 

The presence of inhomogeneous phases in the QCD phase diagram is analyzed within 
chiral quark models that include nonlocal interactions. We work at the mean field level, 
assuming that the spatial dependence of scalar and pseudo-scalar condensates is given by 
a dual chiral density wave. Phase diagrams for Gaussian nonlocal form factors are studied 
in detail and compared with those obtained within the Nambu—Jona-Lasinio model and 
quark-meson approaches. 


I. INTRODUCTION 


Due to the well-known sign problem, present lattice QCD analyses are still not able to pro¬ 
vide fully faithful predictions for QCD thermodynamics at low temperatures and relatively high 
chemical potentials, including the region where the critical point is expected to appear. Thus, our 
knowledge of the strongly interacting matter phase diagram largely relies on the study of effective 
models, which offer the possibility to get predictions of the transition features at regions that are 
not accessible through lattice techniques. In this context, in the last years some works have con¬ 
sidered that the chiral symmetry restoration at low temperatures could be driven by the formation 
of nonuniform phases [l(. One particularly interesting result suggests that the expected critical 
endpoint of the first order chiral phase transition might be replaced by a so-called Lifshitz point 
(LP), where two homogeneous phases and one inhomogeneous phase meet jsj]. This result has been 


obtained in the chiral limit, where the end point 


jecomes a tricritical point (TCP), in the frame¬ 


work of the Nambu—Jona-Lasinio model (NJL) { 3 ]. As it is well known, in this model quark fields 
interact through a local chirally invariant four-fermion coupling. More recently, this issue has also 
been addressed in the context of a quark-meson (QM) model with vacuum fluctuations jj|, where it 
is found that the LP might coincide or not with the TCP depending on the model parametrization. 


In a previous work 


5|] we have analyzed the relation between the positions of the TCP and LP 
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in the framework of nonlocal chiral quark models using a generalized Ginzburg-Landau approach. 
It should be mentioned that nonlocal models can be viewed as extensions of the NJL model that 
intend to represent a step towards a more realistic modelling of QCD. In fact, nonlocality arises 
naturally in the context of successful approaches to low-energy quark dynamics |f], 7], and it has 
been shown [8y that nonlocal models can lead to a momentum dependence in the quark propagator 
that is consistent with lattice QCD results 9-111]. Another advantage of these models is that 
the effective interaction is finite to all orders in the loop expansion, and therefore there is not 
need to introduce extra cutoffs jlij ]. Moreover, in this framework it is possible to obtain an 
adequate description of the properties of strongly interacting particles at both zero and finite 

13-26], The results of Ref. [f| indicate that for all phenomenologically 


temperature/density 


acceptable parametrizations considered the TCP is located at a higher temperature and a lower 
chemical potential in comparison with the LP. Consequently, these models seem to favor a scenario 
in which the onset of the first order transition between homogeneous phases is not covered by an 
inhomogeneous, energetically favored phase. The aim of the present work is to further investigate 
the consequences of the possible existence of inhomogeneous condensates on the thermodynamics 
of nonlocal models by explicitly constructing the associated phase diagrams in the mean field 
approximation. In principle, a full analysis would require to consider general spatial dependent 
condensates, looking for the configurations that minimize the mean field thermodynamic potential 
at each value of the temperature and chemical potential. Since for an arbitrary 3-dimensional 
configuration this turns out to be a very difficult task, even in the case of local models it is 
customary to consider one-dimensional modulations, expecting that the qualitative features of 
the inhomogeneous phases will not be significantly affected by the specific form of the spatial 
dependence carried by the condensates [if]. Due to the additional difficulties introduced by the 
presence of nonlocal quark-quark interactions, here we will consider a simple one-dimensional 
configuration, namely the so-called dual chiral density wave (DCDW) j^, for which the spatial 
dependence of the quark condensates is given by 


(q(x)q(x)) a cos (Q ■ x) 


(q(x)i'y 5 q(x)) a sin(Q • x) , 


(1) 


for both q = u and d quark flavors. Regarding the nonlocal interactions, we will consider the case 
of covariant and instantaneous nonlocal form factors with a Gaussian momentum dependence. 

The article is organized as follows. In Sect. II we present the general theoretical framework 
and propose an ansatz for the bosonic mean fields that leads to the required spatial dependence 
of chiral condensates. The model parametrization is also briefly introduced. Then in Sect. Ill we 
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show the phase diagrams for various parametrizations and discuss the features of the corresponding 
phase transitions. Finally, in Sect. IV we state our conclusions. 


II. THEORETICAL FRAMEWORK 

A. Inhomogeneous condensates in nonlocal chiral quark models 

Let us consider a two-flavor model that includes a four-point coupling between non 
antiquark currents. The corresponding effective action in Euclidean space is given by 


/■ 


G 

Se = I d A x ip(x) (-ifi + m c ) ip{x ) - — j a {x) j a {x) 


ocal quark- 


( 2 ) 


where i/’ is the fermion doublet ^ = (u, d) T and m c stands for the current quark mass in the isospin 
limit. The nonlocal currents j a {x ) are given by 


j a (x) = I d^z Q(z) i/} (x + T a ip (x - ^ 


( 3 ) 


where we have defined T a = (To,T) = (I,ry 5 r), while G(z) is a form factor that characterizes the 
effective interaction. 

The model can be bosonized through the introduction of bosonic fields <h a (x) associated to the 
quark bilinears in Eq. (J3J) 12fl. A standard procedure leads to the Euclidean action 


where 


S E = J d A x d A x' ip{x') D 1 {x',x) tp(x) + ^ J d 4 x <f> a (x) <f> a (x) , 

D~ 1 (x l , x) = 5^\x' - x)(-ip x + m c ) + Q{x' - x) T a <f> a + X ^ 


( 4 ) 


( 5 ) 


We will work within the mean field approximation (MFA), in which the bosonic fields are expanded 
around a real classical configuration d> a (x). After integrating the fermion degrees of freedom one 
obtains the bosonized action 

1 - 


-,(bos) 


/ rf 4 ! / <iV 


Tr log D M l(x'i x ) + 4 > a(®0^a( x ) ^ (4) ( x ' - x ) 


( 6 ) 


where the trace acts on Dirac, flavor and color spaces. 

Let us consider a system in equilibrium at finite temperature T and chemical potential fj,, where 
inhomogeneous phases could be favored. At the mean field level the grand canonical thermodynamic 
potential per unit volume is given by 


WmF — ~ y log Ajf , 


( 7 ) 
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where Z MF is the mean field partition function that arises from the effective action in Eq. Q. If 
the ground state is in general not homogeneous, the quark condensate at a given position x can be 
calculated by introducing an auxiliary static field ip(x). One has 


ip(x)} = - 


5 log Z[(p\ 


5<p(x) 


( 8 ) 


ip=0 


where Z[<p] is obtained from Z MF by changing F (x',x) —> D^ F (x',x) + S^(x' — x) <p(x) in 
the inverse propagator given in Eq. (J5|) , taken at mean field. Moreover, since now parity is not 
necessarily an exact symmetry of the vacuum, one can get in general a nonzero value for the 
condensate {^}{x)Y^{x)) = (^(x) * 75 X 3 '!/’(x)). The latter can be obtained from the partition 
function by adding a term m 4 )(x / — x)i 75 T 3 </?(x) to the inverse propagator in Eq. (J&J) at mean field. 

The thermodynamics can be worked out using the Matsubara formalism. Thus, it is convenient 
to consider the inhomogeneous mean field propagator in momentum space. One has 

-i 


D MF (p',p) = 


(~p + m c ) (27r) 4 5 (4) O'- p) + g 


p + p 


ra $a(p'-p) 


(9) 


where g(j>) and & a (p) are the Fourier transforms of the form factor Q{x) and the mean fields 4> a (x), 
respectively. Since energy is conserved for static mean field configurations, it is also useful to define 
a reduced effective propagator D(p',p ) through 


DmAp'iP) = ( 27 0<KP4 -Pi)D{p',p) . 


( 10 ) 


With these definitions the condensates are found to be given by 

W(J)r.« = -T £ / (0E s TV[UD(p',p)], o = 0,3, (11) 

where the traces are taken over Dirac, flavor and color spaces. Here the fourth component of p' and 
p in D(p',p) is given by p' A = yq = ui n — ip, where p is the chemical potential and ui n = (2 n + 1)ttT 
are the fermionic Matsubara frequencies. 


B. Dual chiral density wave 

We address here the relatively simple situation in which the vacuum is modulated by a dual 
chiral density wave. In this configuration the chiral condensate rotates along the chiral circle, 
carrying a constant three-momentum Q [see Eq. ([T])]. For simplicity, in the following we will 
consider the case of vanishing current quark masses, m c = 0. In this limit, the desired behavior 
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of the chiral condensates can 


ae obtained by considering the following ansatz for the mean field 


configuration & a (p' — P',Q) [281 ]: 

T a $ a {p' - p-Q) = (2vr) 4 (5(p4 -p A ) <$> ^ 1 + s ' 75T 3 ^(3) (p / _ f + s Q) . ( 12 ) 

s=± 

From this ansatz it is evident that the effective propagator will be block diagonal in flavor space, 
thus it can be written as a direct sum of D u and Dd propagators. By calculating the inverse in 
Eq. © we get for the u quark 

; Q) 6® (f -p-Q) A-(p ; Q) 5® (p' - p) 

A + (p ; Q) 6 W(p'-p) 8^\p'-p + Q) 

where A±(p ; Q) and ; Q ) are 2x2 matrices in Dirac space. These are given by 

1 


Du(p',p;Q) = 


(13) 


A±(p-,Q) = 
B± (t ; Q) = 


A (p;Q) 
4>g(t) 
A (t;Q) 


{\pl + ip ± Q) 2 + <t> 2 g(p) 2 ]{ip^ ±p ■ f) + (f> 2 g{p) 2 Q ■ t} 
| [t 2 - Q 2 / 4 + 4> 2 g{t) 2 } 1 - i(tx Q ± t±Q) ■ f] j , 


2 


where 

A (p;Q) = [p 2 -Q 2 /± + ct> 2 g(p) 2 \~ +p 2 Q 2 -{p-Qf . 

The mean field propagator for the d quark is obtained from the previous expressions by 

Dd{p',p;Q ) = D u (p',p;-Q ) . 

In this way, from Eq. (ED we obtain 

(u(x)u(x)} = (d(x)d(x)) = F(Q 2 ) cos (Q ■ x) , 

(u(x)i'j5u(x)) = — (d(x)ij5d(x)) = F(Q 2 ) sin(Q • x) , 

where 

F (q 2 ) = _ 4iV T ^ f d3 p (f>g(p)\p 2 -Q 2 /^ + (t> 2 g(p) 2 ] 


(14) 


(15) 


(16) 


(17) 


(18) 


__( 27r ) 3 A (p-Q) 

with P 4 = uj n — ip. 

If the ground state is assumed to be homogeneous, the mean fields <f> a (x) are uniform. Then, 
from parity invariance one has & a (p) = (27t) 4 <5( 4 ) (p) 5 a o 4>, and the operator in Eq. © can be 
trivially inverted. The corresponding expressions for the condensates are obtained in this case by 
setting Q = 0 in Eqs. (THT) and (fT8l) . namely Q, Q 

/- \ _ /J,\ _ , A T rr ST' f d 3 p g{p) (j) 

(uu) (dd) 4N C T (2tt)3 p 2 + g{pY ^ ’ 

(m7 5 r 3 u) = (<% 5 T 3 d) = 0 . (19) 
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Let us now evaluate the thermodynamic potential in Eq. m for the case of the dual chiral 
density wave. From the mean field partition function Z MF (T, p), using the Matsubara formalism 
we obtain the grand canonical thermodynamic potential per unit volume 

°° f d 3 v 6 2 

w mp (T, p) = -2 N C T J2 /—^ logA (p-Q) +^, (20) 

n =—00 ' 

where once again the fourth component of p in A {p ; Q ) is p^ = ui n ~ ip- Here the integral over p 
is divergent for large momenta. A standard way of regularizing the thermodynamic potential is by 
subtracting the free contribution cuf ree = cu M f(<^ = 0) and adding it in a regularized form. In this 
way one ends up with 

4- — 4- , , re S 
+ 2g + "free i 

where 



= -2 N C T 


OO 

£ 


d 3 p 

(2-/r) 3 


log 


1 + 4> 2 g(p) 


2 0 2 g(p) 2 + 2 (p 2 -Q 2 /4) 
( V 2 + Q 2 / 4) 2 - (p- Q) 2 


LO 


reg 

free 


- N c 


~7tt 2 T 4 
90 


T 2 p? 
3 



( 22 ) 


The mean field values (j) and Q = \Q\ can be obtained by looking for the minimum of w M f through 
the coupled equations 


SUC « F S 

d(f) dQ 


(23) 


A region in which the absolute minimum is reached for a nonzero Q will correspond to an inho¬ 
mogeneous phase. As expected, if chiral symmetry is not dynamically broken (i.e. (f> = 0) the 
regularized thermodynamic potential reduces to the free contribution which does not depend 
on Q. 


C. Model parameterization 

In the chiral limit the model has only one coupling parameter, namely the constant G. In addi¬ 
tion, one has to specify the functional form of the form factor g(p), which requires the introduction 
of some momentum scale A in order to satisfy Lorentz invariance. For definiteness we will consider 
here a Gaussian behavior 


dip) = exp (~p 2 /A 2 ) , 


( 24 ) 


which guarantees a fast ultraviolet convergence of loop integrals. 
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Given the form factor shape, one can fix the model parameters G and A so as to reproduce 
the phenomenological values of the pion decay constant f n and the chiral quark condensate (qq). 
In fact, since we are working in the chiral limit, it is obvious from dimensional analysis that any 
dimensionless quantity turns out to be just a function of the dimensionless combination G = GA 2 , 
while a dimensionful quantity can be written as a function of G times some power of a dimensionful 
parameter, say e.g. f n . According to the recent analysis in Ref. 291, we will take here /^ h = 86 MeV 


and ( qq) ch = —(270 MeV) 3 (superindices stress that the values correspond to the chiral limit), thus 
the “physical” value of G will be that leading to a ratio (— (qq ) 0 ^) 1 ' 3 / f^ 1 — 3.14. In order to check 
the parameter dependence of our results we will consider values for this ratio in the range 2.8 to 
3.5. For = 86 MeV, this corresponds to a shift of at most ~ 30 MeV around the central value 
{-(qq ) ch ) 1/3 = 270 MeV. 

In Fig.[l]the solid lines indicate our numerical results for the parameters G and A that correspond 
to the above range of the ratio (—(< 7 g , ) ch ) 1 / 3 //£ h . In the left panel we show the values of the 
dimensionless parameter G as a function of the ratio (— {qq)^) 1 ^ 3 //^ h , while in the right panel we 
quote the effective cutoff scale A as a function of the quark condensate for the phenomenologically 
preferred value = 86 MeV. The cutoff values are found to be of order ~ 1 GeV, in agreement 
with phenomenological expectations. 

Finally, it is interesting to consider the case of the so-called “instantaneous” form factors, 


which just depend on the three-momentum p 


13 - 


In this case Lorentz symmetry is broken, and 


a spatial cutoff is needed [notice that, in particular, the usual “local” NJL model is obtained by 
setting g(p) = 0(A^ JL — p 2 )\. If we consider once again a Gaussian shape for the form factor, 
namely g{p) = exp(p 2 /A 2 ), and same phenomenological requirements as in the covariant case, the 
corresponding numerical values for G and A are those shown by the dashed lines in Fig. 1. As it 
is discussed in Refs. [3, 3; instantaneous models with soft cutoff functions lead to relative large 
values of the quark condensate, therefore for these models low values of the cutoff are typically 
required. We have taken A = 600 MeV as a lower bound, which implies (— (qq)^) 1 / 3 //^ h > 3.21. 


III. PHASE DIAGRAMS 

Let us consider a hadronic system at finite temperature T and chemical potential p. We will 
discuss the features of the phase transitions in the T — p plane for the nonlocal chiral quark models 
discussed in the preceding section. As stated in the Introduction, in general one can find phases in 
which the chiral symmetry is either broken or restored, and regions in which either homogeneous or 
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Figure 1: Left: values of the dimensionless parameter G for a given ratio (— (qq} 0 ^ 1 ) 1 ^ 3 //£ h . Right: effective 
cutoff scale A for a given value of the quark condensate, with fixed to the phenomenological value (in the 
chiral limit) 86 MeV. Solid (dashed) lines correspond to covariant (instantaneous) Gaussian form factors. 

inhomogeneous phases are preferred. In addition, a general analysis based on the Ginzburg-Landau 
expansion shows that nonlocal models allow the presence of both a tricritical point (TCP) and a 
Lifshitz point (LP), located in different positions Q|. 

We start by discussing the case of Gaussian covariant nonlocal form factors. Our numerical 
results for the corresponding phase diagrams are displayed in Fig. 2, where we show different 
scenarios that may arise if the model parameters lie within the range discussed in the previous 
section. We have chosen four parameter sets, denoted as PI, PII, PHI and PIV, which correspond 
to a pion decay constant = 86 MeV and quark condensate values (—(^) ch ) 1 / 3 = 240, 247, 270 
and 300 MeV, respectively, at zero T and //. The different regions of the phase diagram, as well as 
the corresponding transition curves and critical points, are shown in the left panels of Fig. 2. It is 
seen that in all cases at low temperatures and chemical potentials one finds the usual homogeneous, 
chirally broken (HCB) phase, while for low T and high /r the system lies in an inhomogeneous (IH) 
phase (in the case of PIV, which corresponds to (qq) ch = —(300 MeV) 3 , the onset of the IH phase 
at T = 0 occurs at a chemical potential of about 630 MeV, therefore it is not shown in the figure). 

Let us start by analyzing the case of the parameterization PI (upper left panel of Fig. 2, 
corresponding to a relatively low chiral quark condensate ( qq) ch = —(240 MeV) 3 , and a relatively 
high value of the coupling G ). The different phases are indicated by the shaded areas, while solid 
and dashed lines correspond to first and second order transitions, respectively. At a temperature 
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Figure 2: Left: T—fi phase diagrams for different model parameterizations. Solid (dashed) lines indicate first 
(second) order phase transitions. The dotted line is the lower spinodal corresponding to the homogeneous 
chiral restoration transition, while the dashed-dotted line is a boundary of a region in which there exists 
a local inhomogeneous minimum of the thermodynamical potential. TCP, LP and 3P stand for tricritical, 
Lifshitz and triple points. Right: values of </> and Q as functions of the chemical potential, for T = 0. 


of about 100 MeV and low chemical potentials, it is seen that the system lies in the HCB phase, in 
which chiral symmetry is spontaneously broken. As usual, by increasing ji one finds a second order 
phase transition to an homogeneous phase in which chiral symmetry is restored (HCR phase). If the 
temperature is lowered, the corresponding second order transition curve ends at a tricritical point, 
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beyond which it becomes a first order transition line. Now, by following this line, at a temperature 
X 3 P ~ 20 MeV one arrives at a triple point. For T < T$p, at a given critical chemical potential 
Hc(T) the system undergoes a first order transition from the HCB phase into an IH phase, in which 
chiral symmetry is found to be only approximately restored. On the other hand, if one starts with 
a system in the IH phase and increases the temperature at constant chemical potential, at some 
critical value of T one arrives at a second order phase transition into the HCR phase. As it is shown 
in the figure, the corresponding second order transition line continues beyond the triple point with 
a dashed-dotted line inside the HCB area. The latter represents a boundary of a region in which 
the thermodynamic potential has a local minimum that corresponds to an (unstable) IH phase. 
Finally, in the phase diagram we also show with a dotted line the lower spinodal corresponding to 
the homogeneous chiral restoration transition. 

The previously described first order transition from the HCB to the IH phase is illustrated 
in Fig. 3: left and right panels show contour plots of the mean field thermodynamic potential 
w mI(0; Q) at zero temperature for // = 260 and // = 280 MeV, respectively, which correspond to 
both sides of the transition point fj, c ( 0) = 274 MeV. The plots clearly show the transition from 
an absolute minimum at <j) ~ 340 MeV, Q = 0, to another one in which cj) reduces to about 
50 MeV, while the chiral condensates get spatial dependencies as those given by Eq. (ED, with 
Q ~ 450 MeV. These features are also shown in the upper right panel of Fig. 2, where we quote 
the curves for cj) and Q at T = 0 as functions of the chemical potential. Notice that on the HCB 
side (left panel of Fig. 3) there also exists a local minimum at (0, Q ) ~ (50 MeV, 400 MeV). 

Below the previously discussed phase diagram we show in Fig. 2 the case of parameterization 
PII, in which the quark condensate at zero T and ^ is slightly larger, namely ( qq) ch = —(247 MeV) 3 . 
It can be seen that in this case the IH phase region gets reduced and splits in two: a small “island” 
of IH phase becomes isolated from the large IH phase “continent” found at high chemical potentials 
(see shaded regions in the figure). Then, for PHI and PIV it is seen that the island disappears, and 
the onset of the continent is pushed up to larger values of the chemical potential. The discontinuity 
of Q at this transition for T = 0 becomes increased, as it is shown in the right panels of Fig. 2. 
As stated above, for PIV (lower panels of Fig. 2) it is found that the inhomogeneous phase occurs 
for values of fi out of the region of the phase diagram displayed in our graphs. For comparison, 
in Table 1 we quote the values of the effective cutoffs, the values of cj) and the chiral condensate 
(— {qq)^) 1 / 3 for zero T and /r, and the critical chemical potentials at T = 0, for parameterizations 
PI to PIV. We denote by //'(0) the onset of the IH phase continent, while /r"(0) stands for the 
chemical potential at the transition between the IH phase island and the HCR phase region. All 


11 



<v t° evJ ] 

-1.920E-4 
-1.970E-4 
-2.020E-4 
-2.070E-4 
-2.120E-4 
-2.170E-4 
-2.220E-4 
-2.270E-4 
-2.320E-4 
-2.370E-4 
-2.420E-4 
-2.470E-4 
-2.520E-4 
-2.570E-4 
-2.620E-4 
-2.670E-4 
-2.720E-4 
-2.770E-4 
-2.820E-4 




IGeV"] 

-2.120E-4 
-2.220E-4 
-2.320E-4 
-2.420E-4 
-2.520E-4 
-2.620E-4 
-2.720E-4 
-2.820E-4 
-2.920E-4 
-3.020E-4 
-3.120E-4 
-3.220E-4 



Q [MeV] 


Q [MeV] 


Figure 3: Contour plots of the thermodynamic potential u%f for a nonlocal chiral quark model at zero 
temperature and finite chemical potential, close to the first order transition between HCB and III phases. 
The plots correspond to parameterization PI, for chemical potentials fj, = 260 MeV (left) and p, = 280 MeV 
(right). 


(- 

-( qq) ch y /3 

A 

<t> 

Mc(0) 

m"(o) /4(o) 

PI 

240 

808 

338 

274 

- 

PII 

247 

863 

315 

266 

288 295 

PHI 

270 

1045 

264 

249 

470 

PIV 

300 

1295 

227 

236 

629 


Table I: Effective cutoff A, chiral condensate and mean field (f> at zero T and /i, and zero-temperature critical 
chemical potentials for parametrizations PI to PIV. All values are given in MeV. 

values are given in MeV. 

It is worth pointing out that, for these models, the would-be Lifshitz point (i.e. the point where 
the HCB and HCR phases would meet the IH one along the second order phase transition lines) 
is hidden inside the HCB phase region: around the would-be LP, the HCB phase turns out to be 
energetically preferred in all cases considered. Instead, as it is indicated in the upper left panels of 
Fig. 2, a triple point can be found in the case of PI and PII. It is also worth mentioning that the 
second order phase transition curves, as well as both the TCP and would-be LP, can be calculated 

nn 

for these models through a quite precise semianalytical approach [51. 1181]. 

The characteristics of the phase diagrams can be compared with those obtained within the NJL 
and the quark-meson model, which have been recently analyzed in this context. As stated in the 
Introduction, in the NJL the TCP and LP are coincident, whereas in the QM this can be the case 
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Q 


or not, depending on the parameterization. In some cases it is shown [J] that the HCB-HCR second 
order phase transition ends at a Lifshitz point, while the TCP appears to be hidden into the IH 


region (i.e., the opposite situation to t 


according to the analyses in Refs 


■ HQ 


lat found in our models). It is interesting to notice that, 


311 ]. for both the NJL and QM models some parameteri- 


zations lead to phase diagrams that show IH “continents” that extend to arbitrarily high chemical 
potentials. In fact, it is a matter of discussion whether the presence of these continents arises just 
as a regularization artifact. We stress that in nonlocal models the ultraviolet convergence of loop 
integrals follows from the behavior of form factors, which effectively embrace the underlying QCD 
Re form factors can be fitted from lattice QCD calculations for the effective 


quark propagators 



re 

I 

26 


261] h The fact that various quark models including different regularization 


procedures lead to similar qualitative features of the phase diagram seems to indicate that these 
features are rather robust. However, it is necessary to mention that we have not considered the 
effects of color superconductivity, which are expected to be important at intermediate and large 
chemical potentials and could have a significant impact in the phase diagram. 

Finally, we address the case of instantaneous form factors mentioned in the previous section. 
For these parametrizations the qualitative features of the phase diagrams are found to be basically 
the same as those discussed above. The main difference is that similar diagrams correspond to 
models leading to larger values of the quark condensates: for a model with (— (qq)^) 1 / 3 = 270 MeV 
one obtains a phase diagram similar to that of parametrization PI for the covariant case, while the 
small “island” of inhomogeneous phase arises when the corresponding condensate is (— (qq)^) 1 / 3 ~ 
285 MeV. For larger absolute values of the quark condensates the island disappears and the onset 
of the inhomogeneous phase is pushed up to larger values of the chemical potential, just as in the 
case of covariant form factors. 


IV. SUMMARY AND CONCLUSIONS 

In this work we have analyzed the possible existence of inhomogeneous phases in the context of a 
simple version of nonlocal SU(2) chiral quark models in the chiral limit. For simplicity, only the one- 
dimensional modulation associated to a dual chiral density wave (DCDW) has been considered. 
In this framework, different parametrizations of the nonlocality, including both covariant and 
instantaneous form factors, have been investigated. 

For all studied scenarios it is seen that the sizes of inhomogeneous phase regions show a rather 
strong dependence on model parameters. In all cases we find the existence of a tricritical point, 
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while, keeping f n fixed, for high values of the dimensionless coupling G (low absolute values of the 
chiral condensate) we find at low temperatures a first order transition between the homogeneous 
chirally broken phase and the inhomogeneous phase. These phases and the homogeneous chirally 
restored one meet then at a triple point. On the other hand, for lower values of G the onset of the 
inhomogeneous phase is pushed up to higher chemical potentials and the triple point disappears. 


As in previous analyses made in the framework of NJL and QM models 


311 ]. the inhomoge¬ 


neous “continents” in our phase diagrams extend to arbitrarily high chemical potentials. Thus, 
their existence seems to be a rather robust prediction of this type of quark models. It should be 
mentioned, however, that effects of color superconductivity, which are expected to be important at 
intermediate and large chemical potentials, have not been included in these works. In this sense, it 
is clear that to clarify this issue color superconducting interaction channels have to be incorporated 
in future calculations. In addition, the role of vector channels within the framework of the non¬ 
local models deserves further investigation as well. Finally, it would be interesting to explore the 
possibility of going beyond the DCDW ansatz used in the present nonlocal scheme by considering 
more general one dimensional modulations. 
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